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Abstract. The variability and rotation of ultra cool dwarfs (UCDs) provide important information on the atmospheres and 
evolution of these very low mass stars and brown dwarfs. As part of an ongoing program to investigate this, the projected 
rotation periods, vsini, derived from high resolution VLT/UVES spectroscopy via cross correlation are presented for 16 field 
UCDs (M9V-L7.5V). This doubles the number of L dwarfs for which vsini has been measured. All targets are found to 
have vsin; between 10 and 40km/s confirming that L dwarfs are rapid rotators. Radial velocities have also been measured 
to a precision of 1-2 km/s. From the random distribution of the rotation axes, i, and theoretically predicted radii, one-sided 
confidence intervals are placed on the rotation periods of individual objects. These are compared with published period data 
obtained from photometric monitoring programs. From this, the period of 31 hrs for the L0.5 dwarf 2M0746+2000 published 
by Gelino et al. 1 2002 1 may be ruled out as the rotation period. The period of 11.2 ± 0.8hrs for the LI. 5 dwarf 2M1 145+2317 
obtained by Bailer- Jones & Mundt 12001 1 is consistent with the present vsini results so is plausibly the true rotation period. 
' The inclination of the rotation axis is constrained to be i=62°-90° with an expectation value of 76°. Alternatively the data set a 

lower limit on the radius of 0. 1R Q , which is within the range of radii predicted by models for brown dwarfs older than 0.5 Gyr. 
Similarly, the period of 2.7 ±0.1 hrs detected by the same authors for 2M1334+1940 is also confirmed as the likely rotation 
period; the inclination is i=27 c -44°(< ( >= 34°). Where no variability or period was detected by the monitoring programs the 
likely reason is low contrast modulating surface features. However, in three cases variability but no period was detected, even 
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though the likely rotation period range inferred from v sin i lies within the timescale to which the monitoring was sensitive. This 
reinforces the 'masking hypothesis' of Bailer- Jones & Mundt (2001 1, the idea that the evolution of photospheric features on 
timescales shorter than the rotation period obscure the regular modulation of the light curve. As has been previously discussed, 



a likely candidate for such features is inhomogeneous dust clouds. 
Key words. 

1 . Introduction I1999I I. In particular, global dust structures ('clouds') may play 

a significant role. 

Understanding the nature of low mass stars and brown dwarfs 

requires a detailed characterization of their observable atmo- Significant information on UCDs can be obtained from 
spheres. Specifically, the determination of luminosities, chem- ^ observable temporal changes, and several dedicated moni- 
ical compositions and ages ultimately depends on adequate toring programs have now uncovered good evidence for optical 
knowledge of the radiative and convective transport mecha- and infrared variability (Bailer-Jones & Mundt 1999 Terndrup 
nisms. While our understanding of these mechanisms is rea- et al.H222 Tinney & Tolleyll999] Nakajima et al-ESDI Bailer- 
sonable (if incomplete) for some types of stars, it is compara- Jones & Mundt EMU Martin et al. [2001 Bailer-Jones ES2 
lively poor for ultra cool dwarfs (UCDs), late type M, L and T Burgasser et al. 2002 Clarke et al.|2002ajj2002b GelinoEE2 



dwarfs. This is not least because at the low effective tempera- Gelin o et al - EQQ1 Bailer -Jones & Lamm 2003 Clark e et al 

tures involved (<3000K), solid dust particles form in signifi- EQQ3 Enoch et al. |20 j)3| Jo ergens et al. |2003] Koen |2003j 

cant numbers and varieties, and these have a major impact on Zapatero Osorio et al| 2003J and Ha variability (Hall |2002J 

the structure of the atmosphere and hence on the observable Mochnacki et al. |2002J Liebert et al. |2003). Of the 80 or so 

properties of the star (e.g. Allard et al.ESD Burrows & Sharp UCDs monitored by these groups, about 30 show evidence 

for photometric variability. While this figure depends on what 

* Based on observations obtained at the VLT, European Southern one takes as sufficient evidence for variability, at least half are 

Observatory, Chile, program 65.L-0199 convincingly variable on the scale of a few tens of millimag- 

** Emmy Noether Fellow of the Deutsche Forschungsgemeinschaft nitudes. There are indications that this photometric variabil- 
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ity is not caused by a simple rotational modulation of a non- 
uniform photosphere. One candidate is an inhomogeneous dust 
layer evolving under the influence of convection on time scales 
shorter than the rotation period (see the introduction to Bailer- 
Jones 2002 for a discussion). Understanding this phenomenon 
is important, not only for what it tells us about atmospheric pro- 
cess under these interesting physical conditions, but also for the 
implications it has concerning conclusions drawn from single 
epoch observations of UCDs. 

In this article I address the specific issue of how the pho- 
tometric variability and period data relate to UCD rotation, 
as measured by line broadening from high resolution spec- 
troscopy. This provides the projected rotation speed, v sin i, so 
using simple statistical arguments and theoretical predictions 
for UCD radii, constraints can be placed on likely rotation pe- 
riods. These (plus other published vsim data) are compared 
with the variability data to confirm or refute possible rotation 
periods and investigate the 'masking hypothesis'. 

Previously published work has established that UCDs are 
fast rotators: of the 17 hitherto observed L dwarfs, 16 have 
vsin/ in the range 10^-0 km/s and one (Kelu-1) has 60km/s 
(Basri et al. |2000 Mohanty & Basri 120031 Schweitzer et al. 
2001 1. In the present work, v sin i values are determined for 16 
UCDs (15 L dwarfs and 1 M dwarf), 14 of which have no previ- 
ously published measurement. Moreover, conservative ranges 
for v sin i are established which reflect the dominant uncertain- 
ties in measuring v sin i by cross correlation. For the majority of 
UCDs, the complete v sin i range lies between 10 and 30km/s, 
and there is no measurement below about lOkm/s, thus con- 
firming the rapid rotation of UCDs. For a discussion of UCD 
rotation and its relation to chromospheric activity and possible 
dynamo mechanisms, the reader is referred to Mohanty & Basri 
(2003. 

2. Data acquisition and reduction 

A target list of UCDs ranging in spectral type from M9 to L8 
was assembled from the published UCDs available at the time 
of the observations (April 2000). They were selected preferen- 
tially for brightness and for good observability from the VLT 
and to cover a range of spectral types. They were not selected 
based on any known v sin i values (none were available at the 
time), although some L dwarfs monitored by Bailer-Jones & 
Mundt \ 19991120011 were specifically included. 

2.1. Instrumentation and observations 

High resolution spectroscopy was obtained with the UVES 
echelle spectrograph mounted on the UT2 8.2m VLT telescope 
at Cerro Paranal, Chile, on 26-28 April 2000. The red arm of 
this instrument was used with the CD4-prot cross dispersion 
grating and a slit width of 1". This provides the wavelength 
range 6440-10250 A (orders 94-60) at a resolution of 38700, 
yielding a FWHM of 7.8 km/s. The detector is a mosaic of two 
CCDs with pixels binned to provide a sampling of 2.4km/s/pix 
in the dispersion direction and 0.35"/pix in the spatial direc- 
tion. The slit length was 10" and was oriented vertically (i.e. 
parallel to the direction of maximum atmospheric dispersion). 



Over the course of two and a half nights, spectra were ob- 
tained of the 16 UCDs listed in Table|2] Exposure times ranged 
from 1x15 min for the brightest targets to 3x50 min for the 
faintest targets. The seeing was typically 0.8" (0.5"— 1.0") at 
airmasses of up to 2.0. Spectra were also obtained several times 
with the same instrument settings for v sin i and v ra d standards, 
namely the bright M dwarfs G1402, G1406 and G1876. The 
spectrograph is mounted on the Nasmyth platform so under- 
goes limited movement during the course of the observations. 
The data reduction procedure verified the stability of the instru- 
ment. While wavelength calibration spectra and flat fields were 
taken at regular intervals each night, it was found that a single 
set of calibration frames for each night was sufficient. 

2.2. Data reduction and spectral extraction 

The objective of the data reduction is to produce wavelength 
calibrated one-dimensional spectra for each echelle order, free 
of relevant instrumental and telluric signatures as preparation 
for cross correlation. The basic data reduction was carried out 
using the IRAF 1 package, and consists of the removal of scat- 
tered light, flat fielding, spectral extraction and wavelength cal- 
ibration. 

The two-dimensional data format (i.e. cross dispersed spec- 
tra) will be referred to as an "image". Scattered light in the 
spectrograph contributes a two-dimensional additive pattern to 
the dispersed light and generally needs to be removed from 
both science frames and flats. This was done by tracing the 
orders and defining the inter-order regions. A two-dimensional 
fit is made to these regions and this fit subtracted from all im- 
ages. As dark frames, biases and overscan regions showed no 
features, a further explicit subtraction of these was not neces- 
sary. 

Flat fielding is performed to remove small scale detector 
variations and fringing. 2 A sequence of individual flat field im- 
ages is combined. This master flat was then rectified using the 
apflatten task, which makes fits to each order to remove the 
pseudo-continuum and grating blaze function. The result is a 
two-dimensional flat field preserving only high frequency vari- 
ations which are then removed from the science frames by divi- 
sion. The global variations (pseudo-continuum, blaze function 
etc.) are removed later immediately prior to cross correlation. 

The spectra are extracted and the sky subtracted using the 
apsum package. A one-dimensional optimal extraction was per- 
formed with pixel cleaning based on expected noise statistics. 

1 the Image Reduction and Analysis Facility, provided by the 
National Optical Astronomy Observatories (NOAO). 

2 In broad band direct imaging, fringing (monochromatic interfer- 
ence in the CCD) is generated by night sky line emission sources. As 
discussed in Bailer-Iones & Mundt 1 2001 1 this produces an additive 
pattern over the whole CCD which must be subtracted to perform cor- 
rect photometry (and not divided, as is often thought, because it does 
not modulate the star light). With high resolution spectroscopy the sit- 
uation is different. Here the grating itself creates an independent nar- 
row band source at every point along the dispersion direction which in 
turn generates fringing at each point in the CCD. All incident sources 
therefore create the same fringe pattern at a given wavelength, so the 
fringes must be divided out using a flat field. 
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Because many of the L dwarfs are faint, extraction apertures 
are not traced: they are transferred from the order definition 
image and simply recentered and resized to accommodate both 
the position of the sources in the slit and the seeing. 

A wavelength calibration was established using the internal 
ThAr lamp. Lines were iteratively identified and a fit produced 
using the ecreidentify package. The fit was a two-dimensional 
Chebyshev polynomial based on the spectral order numbers 
(using a third order polynomial along the dispersion axis and 
second order one perpendicular to it). The RMS of the fit was 
0.01 pixels, or 25m/s. This was applied to the science spectra 
using the dispcor package. Changes in temperature or pressure 
or seismic disturbances can alter the zero point of this calibra- 
tion during the course of a night. The calibration was there- 
fore checked in all science images against the positions of the 
OH emission lines from the Earth's atmosphere. These demon- 
strated that the zero point shift of the instrumental calibration 
was below the RMS of the fit. Of course, additional zero point 
shifts occur in the applied wavelength scale if a star is not per- 
fectly centered in slit, and could be as large as the projected slit 
size. In practice, though, autoguiding keeps the star centered in 
the slit to within about a tenth of the FWHM, giving an upper 
limit on the systematic zero point shift (and hence systematic 
radial velocity error) of about 1 km/s. This does not effect the 
rotational velocities. 

Candidate wavelength regions for the cross correlation 
were selected on the grounds of being relatively free of tel- 
luric absorption, as identified via inspection of stellar spectra 
with few intrinsic features (a rapidly rotating B9V, a B5V and 
an AOIII) acquired for this purpose. As a result, just those or- 
ders (or parts of orders) listed in Table ^ are retained for all 
targets and templates. Regions containing strong telluric emis- 
sion lines (which, on account of noise, are only poorly removed 
by sky subtraction) were likewise removed. While performing 
the cross correlations on individual objects, all orders were vi- 
sually inspected and any obvious emission lines, cosmic rays 
or bad pixels were masked. Not all orders in Table[0were used 
to form the final values of vsin/ or v ra d: orders were selected 
for each star separately based on the quality of the cross corre- 
lation function (see below). Example spectra for three objects 
are shown in Fig.^ The signal-to-noise ratio (SNR) varies for 
different orders and different targets, but in order 70 it is 8-25 
per extracted pixel (0.05 A; Cf. FWHM of a resolution element 
of 0.23 A in order 70). Orders with SNR lower than a few were 
disregarded; for this reason the bluest orders were often not 
used. 

3. Radial velocities 

Radial velocities are determined via cross correlation of 
the UCD spectra with the spectrum of an observed star of 
known radial velocity. This was done in two stages. The ini- 
tial template used was G1406, an M4V star with v ra d re- 
ported as 19.48 ±0.03 km/s by Nidever et al. (|2002i and 
19.18 + 0. 1 1 km/s by White & Basri ( 120031 . 1 adopt the former. 
The cross correlation was performed using the IRAF package 
fxcor, with the centre of the cross correlation peak calculated 
using the 'centerld' option. The pseudo-continuum of each or- 



Table 1. Candidate orders or partial orders retained for the 
cross correlation. In individual cases, some orders produce 
poor correlation functions in either the cross correlation of the 
template against the target or, for v sin i, the template against 
the spun up template (calibration), in which case the order is 
not used. Orders 90-94 were not used with the L dwarf rota- 
tion template, and only orders 60-72 were use in determining 
v ra d, m both cases for SNR reasons. 



order 


^min./ A 


'I max./ A 


60 


10095.0 


10250.0 


61 


9931.0 


10080.0 


62 


9850.0 


9917.0 


69 


8786.0 


8903.0 


70 


8661.0 


8775.0 


71 


8540.0 


8651.0 


72 


8419.0 


8530.0 


75 


8087.0 


8130.0 


76 


7979.5 


8078.0 


77 


7877.0 


7973.0 


78 


7777.0 


7869.0 


79 


7699.0 


7770.0 


81 


7490.0 


7577.0 


82 


7399.0 


7484.0 


85 


7139.0 


7165.0 


86 


7054.0 


7133.0 


89 


6820.0 


6866.0 


90 


6745.0 


6815.0 


91 


6670.0 


6739.0 


92 


6598.0 


6665.0 


94 


6460.0 


6523.0 




B660 B680 8700 8720 B740 8760 8780 

wavelength / Angstroms 



Fig. 1. Spectra (order 70) of an M dwarf v sin i and v ra d stan- 
dard (top), an L1.0 dwarf (also the L dwarf 'standard', mid- 
dle), and an L3.0 dwarf (bottom). The spectra have been con- 
tinuum subtracted, scaled to a common vertical range and offset 
on the vertical axis. The exposure times for these single spec- 
tra were 3 min, 15 min and 40 min respectively, which achieve 
mean SNR per resolution element around the centre of the or- 
der of 150, 15 and 1 1 respectively. 
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der is rectified (fit and divided out) using a fifth order cubic 
spline. 

The cross correlation function was sometime quite noisy, 
in particular for the bluer spectral orders and for late L spectral 
types, resulting in ambiguous peak centering or even identifica- 
tion. For this reason, the bright L1.0 dwarf 2M1439+1929 was 
adopted as a 'secondary' radial velocity standard, and only the 
orders in the 'red' CCD (orders 60, 61, 62, 69, 70, 71 and 72 in 
TableQ} are used. The radial velocity of 2M1439+1929 is de- 
termined to be -26.3 ± 0.50 km/s from cross correlation against 
G1406. (The bright L0.5 dwarf 2M0746+2000 could also have 
been used for this purpose, but was not on account of its bi- 
narity; section^}. The radial velocities for all other UCDs are 
determined by cross correlation against 2M1439+1929 and are 
listed in Table |2 the value given is a weighted mean of the 
results in the different orders. The uncertainty is the standard 
deviation in this mean (including the uncertainty in the tem- 
plate velocity), apart from for a few cases in which the peak 
centering uncertainty was larger, in which case this is reported 
instead. 

Three of the results are dubious (marked with a ":") because 
of discrepancies between orders and/or significant centering 
uncertainties due to highly rotationally broadened spectral 
lines. The latter is particularly applicable for DEI 159+0057 
and DEI 43 1-1953. Two of the targets have v lat j determinations 
in the literature. Basri et al. (2000) obtained v ra( j = 29+1 km/s 
for 2M1439+1929, which differs from my result by 2.5cr. Reid 
et al. J2002L on the other hand, found v lad = 54.1 ±0.80 km/s 
for 2M0746+2000, in excellent agreement with my result. 
Recall that my results could have an additional systematic error 
of up to 1 km/s (section l2~2l . 

I should note that Basri et al. |2000| have observed G1406 to 
have a significantly different radial velocity of around 40 km/s 
on one epoch (June 1997), although they found a constant 18- 
19 km/s at several other epochs. They argue that this anomaly 
is real, perhaps due to stellar pulsations or a companion in a 
highly elliptical orbit. I therefore checked the radial velocity of 
G1406 against two other standards observed, G1402 and G1876. 
These give 18.8 ±0.1 km/s and 18.2 ±0.5 km/s respectively 
(adopting v lad = - 1 .04 ± 0.03 km/s for G1402 from Nidever et al. 
120021 and v rad = -2.0 ± 0.5 km/s for G1876 from Delfosse et al. 
I1998> . This is slightly lower than the adopted value for G1406. 
Although the difference is significant compared to the random 
errors, it is at a level where systematic errors could dominate. 
Moreover, it has no statistically significant impact on the de- 
rived radial velocity for the secondary template 2M 1439 +1929: 
cross correlation against the G1402 and G1876 spectra yields - 
26.6 ± 1.3 km/s and -27.3 ± 1.2 km/s respectively, compared to 
the adopted value of -26.3 ± 0.50 km/s. 

4. Rotational broadening via cross correlation 

UCD rotation velocities are determined via a cross correlation 
of the UCD dwarf spectra (the 'targets') against spectra of stars 
with essentially zero rotation velocity (the 'templates'). Under 
the assumption that the line broadening of the UCD is dom- 
inated by the rotational broadening (and that the line broad- 
ening of the template is comparatively small), the width of the 



cross correlation function is a measure of the v sin i of the UCD, 
which may then be calibrated as discussed below. This follows 
the method used by Tinney & Reid ( 1998 1, Mohanty & Basri 
(2003 ) and White & Basri ( 2003 ). The cross correlation is done 
on each order separately and the results combined to obtain the 
final v sin i measure. 

The two primary templates are G1402, an M4V with 
vsin/<2.3km/s, and G1406, an M5.5V with v sin ;'<2.9 km/s 
(Delfosse et al. I1998> . The procedure implicitly assumes that 
the UCDs have similar spectra to the templates, ideally differ- 
ing only in v sin i. This is of course not the case, and the dif- 
ferences between mid M dwarfs and early/mid L dwarfs could 
bias the results. But in the absence of appropriate L dwarfs with 
independently determined rotation velocities this is the best we 
can presently do with the cross correlation method. (An alterna- 
tive method is to examine individual lines according to a phys- 
ical model of line formation and achieve a best fit against a 
grid of model spectra, as done, for example, by Schweitzer et 
al. 2001 , although this was not very sensitive to v sin z). In some 
cases, poor fits to the cross correlation peak occur, and such fits 
are rejected. The impact of spectral type mismatch is assessed 
via a 'bootstrap' method, namely by using the slowest rotating 
L dwarf as a secondary template (see section[5}. 

The cross correlation is done again with fxcor, but inde- 
pendently of the radial velocity determination. The main peak 
of the cross correlation function is fit with a Gaussian and its 
FWHM measured. If a good fit cannot be obtained, or if the 
peak is ambiguous, then the order is not used. 

The main uncertainty in performing this fit is determin- 
ing the 'background' for fitting the Gaussian (i.e. the effec- 
tive zero level in the cross correlation function). For this rea- 
son, three Gaussians were fit corresponding to the best fit, the 
minimum FWHM fit and the maximum FWHM fit. The fits 
were performed interactively, with the maximum and mini- 
mum representing conservative limits. The range between these 
comprises a generous assessment of errors in the fitting proce- 
dure. Although the peaks sometimes deviated from a Gaussian 
shape (especially at larger rotation velocities), a Gaussian was 
nonetheless felt to be the best overall parametrized form among 
those tested. 

These FWHM values are converted to v sin i values through 
the following calibration process. The template spectral orders 
are artificially spun up by convolving their spectra with the 
rotation profile given by Gray J 19921 eqn. 17.12) with vsin/ 
spanning 10 to 100 km/s in steps of 5 km/s. A limb darkening 
parameter of e = 0.6 is assumed in this profile. The spun up 
spectra are cross correlated with the non-rotating template and 
the FWHM of the peak is measured in the same way as de- 
scribed above. (Here the difference between the FWHM of the 
minimum and maximum Gaussian fits is much smaller than in 
the UCD case, so is neglected.) A linear least squares fit is then 
obtained between the FWHM and v sin /, separately for each 
order. It was found that the Gaussian was a poor fit to this cross 
correlation function above vsin/ = 70 km/s, so the fit is only 
made in the range 10-70 km/s. For a few orders a good calibra- 
tion could not be obtained so these were not used. 

When this calibration is applied to determine v sin i for each 
order and these combined (next section), we must pay attention 
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to the minimum value of v sin i to which this procedure is sen- 
sitive. Following Tonry & Davis ( 1979 ), the measured width of 
the cross correlation function is made up of a number of com- 
ponents 

^meas = ^rot + ^nat + 2 ^Lt ( 1 ) 

where <x rot is the rotational width, cr nat is the natural, or intrin- 
sic, line width in the UCD, and <x inst is the instrumental broad- 
ening (two contributions as two spectra form the cross corre- 
lation function). I assume that the rotational and intrinsic con- 
tributions from the template are negligible. Because the terms 
are additive, then with 'perfect' data, any non-zero rotational 
broadening could be detected even with non-zero <x mst . In real- 
ity, however, small broadenings cannot be detected due to noise 
and template/target mismatch. I make the somewhat ad hoc as- 
sumption that we can only detect cr ro t if it exceeds the other 
broadening contributions. If we further assume that the UCD 
natural broadening is negligible (reasonable if resonance lines 
and gravity sensitive lines are avoided), then the minimum de- 
tectable value of cr 10t is equal to V2cr inst . The FWHM of the in- 
strumental broadening is set by the slit width of 7.8 km/s (equal 
to 1"; when the seeing was below this - as was often the case 
- the instrumental broadening is reduced). Note that the full 
width of a rotational profile corresponds to twice the rotational 
velocity (one half of the line is created by the blueshifted ap- 
proaching limb of the star, the other half by the redshifted re- 
ceding limb). Thus the minimum detectable vsin/ is therefore 
V2 x 7.8/2 = 5.5 km/s. Any derived values of v sin i below this 
limit are dropped from the calculations in the next section. If 
this limit has been overestimated, then we would potentially 
drop too many measures, overestimating the combined v sin i, 
(It turns out that I only drop one value for being below this 
limit, so this is not significant.) If this limit has been underes- 
timated, then additional values may need to be dropped, poten- 
tially raising v sin i. 

5. Projected rotational velocities 

The calibrations described in the previous section are applied 
to each order to determine a v sin i for each of the three fits. The 
different orders are then combined to give the mean, minimum 
and maximum values as listed in Table |2 The mean value is 
a weighted mean of the 'best fit' v sin / values described in the 
previous section, using a weight of 1,2,3 or 4 depending on the 
quality of the fit to its cross correlation peak. Obvious outliers 
are clipped when forming this mean (7% of all orders/spectra 
are clipped in this way). The minimum and maximum values 
listed in Table|2]are likewise a weighted mean of the minimum 
and maximum derived v sin i values. 

As discussed in the previous section, the cross correlation 
method implicitly assumes that the template and target spectra 
differ only in the rotation velocity, yet the templates are mid 
M dwarfs whereas the targets range in spectral type from M9 
to L5 with one L7.5. This potentially biases the v sin i determi- 
nations. To partially overcome this, the entire cross correlation 
and calibration procedure is repeated using the slowest rotat- 
ing L dwarf as a template. The relevant star is 2M1439+1929, 
which was also used as a template by Basri et al. (2000 1 for the 



same reason and purpose. Fortunately, this LI .0 dwarf is bright 
enough to use as a template and is a reasonably early L type 
such that its cross correlation against the mid M dwarfs was 
probably not too erroneous. Of course, with v sin z'=l 1.2 km/s, 
this L dwarf is rotating faster than the M dwarf templates. 
Using a rotating template will generally underestimate v sin 2, 
because when we spin up the templates for the calibration, a 
given vsin/ corresponds to a higher FWHM than would be 
the case with a non-rotating template. The derived v sin / val- 
ues using 2M1439+1929 as the template are shown in Table[2] 
Other than in one case (and that is a single measure of lowest 
quality), there is actually no significant tendency for v sin / to 
be lower (or higher) than the values derived with the M dwarf 
templates. Indeed, the agreement between the values is good, 
deviating by less than 2 km/s in 10 out of 14 cases. In two cases 
(DE1431-1953 and DEI 159+0057) the disagreement is more 
than 10 km/s, but in both of these cases one of the measure- 
ments is based on only a single (poor quality) measurement. In 
the remaining two cases (2M0913+1841 and SD1203+0015) 
the discrepancy is around 5 km/s but the minimum to maxi- 
mum ranges just about overlap (and both stars are among the 
faintest observed). Thus although the vsirn' calibration using 
2M1439+1929 has used the M dwarf calibration to identify 
this template, the spectra and fitting procedures are indepen- 
dent. The fact that both templates give consistent v sin i ranges 
indicates that the impact of template mismatch it not major, at 
least not when using multiple orders across a wide wavelength 
range. 3 Fig. ^ further shows that there is relatively little mis- 
match between even an M5.5V and L3.0V. Only three of the 
UCDs are later than this. For the latest type (L7.5V), no ac- 
ceptable cross correlation function could be obtained against 
the M5.5V, perhaps indicative of commencing mismatch prob- 
lems. In the rest of this article I adopt the 2M1439+1929 cross 
correlation values, except of course for 2M1439+1929 itself 
and except for DEI 159+0057 because it has only a single poor 
quality measurement against this template. 

The most striking result from Table |2 is that L dwarfs 
are relatively fast rotators, considerably faster than is typically 
found for M dwarfs. They all have vsin/ in the range 11- 
37 km/s. Moreover, the lowest value of vsin; is 11. 2 km/s, and 
the lowest minimum value is 9.4 km/s, well above the minimum 
detection limit of 5.5 km/s (section^. This agrees with the re- 
sults of Mohanty & Basri ( 2003 1, who found no L dwarf in 
their sample of 13 to have vsin / below lOkm/s. Two specific 
targets, 2M0746+2000 and 2M1439+1929, have vsirn' deter- 
mined by others: The former has 24 ± 5 km/s according to Reid 
et al. (2002 1 and 20 + 10 km/s according to Schweitzer et al. 
The latter is derived to have 10 ± 2.5 km/s by Basri et al. (2000l 
and 10 + 2 km/s by Mohanty & Basri (2003). These agree well 
with the values and ranges derived in the present work. 



3 This consistency could be misleading if use of the L template 
gives a positive bias to the derived vsin / (compared to the M tem- 
plate) which is then offset by a negative bias from having used a rotat- 
ing L template. However, it is not obvious that the later-type template 
should give a positive bias and, moreover, the two biases would have 
to contrive to have the same magnitude. 
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Table 2. v sin i values for ultra cool dwarfs observed in this survey. The first column lists the full name of each target: an 
abbreviated version is used elsewhere in this paper, where 2M=2MASS, DE=DENIS and SD=SDSS. The second and third 
columns list the spectral types and I magnitudes taken from the literature. (I magnitudes given to only one decimal place have 
been estimated.) Columns 4-7 and 8-11 list the values of vsin/ deduced via cross correlation with the M dwarf templates and 
L dwarf template (2M1439+1929) respectively, "mean" (i.e. vsin/) is a weighted mean formed by combining measurements 
from individual spectrograph orders and/or exposures: N orders are retained in the final calculation, "min." and "max." are 
the minimum and maximum values of v sin i similarly averaged over all orders, and so encompass generous estimates of the 
uncertainty in fitting the peak of the cross correlation function. The final column gives the radial velocities relative to the solar 
system barycentre. In addition to the random errors given for these, there may be an additional systematic error from the zero 
point uncertainty (due to possible slit decentering) of no more than 1 km/s. A colon, indicates an uncertain value. 



Name 




SpT 


/ 








v sin i 1 km/s 








Vrad / km/s 










M dwarf template 




2M1439+1929 template 














mean 


min. 


max. 


N 


mean 


min. 


max. 


N 






DENIS-P 


11431-1953 


M9.0 


17.47 


47.4 


44.1 


52.6 


1 


37.1 


33.7 


41.6 


4 


-11.4± 


2.8 


DENIS-P 


Jl 159+0057 


LO.O 


17.32 


35.2 


31.9 


39.7 


3 


74.5 


69.1 


83.5 


1 


-3.0 ± 


1.5 : 


2MASSI 


10746425+200032 


L0.5 


15.11 


27.3 


25.6 


30.6 


6 


25.8 


23.0 


28.8 


9 


54.1 ± 


0.9 


2MASSW 


J1412244+163312 


L0.5 


17.1 


17.3 


14.7 


19.2 


6 


16.4 


13.9 


19.2 


4 


6.9 ± 


0.8 


2MASSW 


J1439284+192915 


L1.0 


16.12 


11.2 


9.6 


12.8 


9 










-26.3 ± 


0.5 


DENIS-P 


J1441-0945 


L1.0 


17.32 


17.4 


14.8 


20.6 


4 


15.9 


13.3 


19.3 


6 


-27.9 ± 


1.2 


2MASSW 


Jl 145572+23 1730 


L1.5 


18.62 


12.5 


10.5 


14.4 


9 


12.7 


10.5 


14.8 


6 


3.7 ± 


0.9 


2MASSW 


J1334062+194034 


L1.5 


18.76 


25.2 


22.1 


28.1 


4 


25.4 


21.5 


29.9 


5 


-4.3 ± 


2.1 : 


2MASSI 


J1029216+162652 


L2.5 


17.9 


28.3 


23.6 


32.5 


7 


28.0 


22.8 


33.1 


8 


-29.2 + 


4.0 


DENIS-P 


J1047-1815 


L2.5 


17.75 


16.3 


15.0 


18.3 


12 


15.0 


12.2 


17.6 


11 


6.0 ± 


0.8 


2MASSW 


J0913032+184150 


L3.0 


19.07 


15.0 


11.9 


17.6 


11 


20.3 


18.0 


22.6 


3 


28.4 + 


2.5 : 


SDSSp 


J120358. 19+001550.33 


L3.0 


18.88 


27.6 


24.1 


32.3 


7 


31.7 


25.6 


36.8 


8 


-2.7 ± 


2.3 


2MASSW 


J1615441+355900 


L3.0 


18.1 


12.1 


10.7 


13.8 


5 


12.8 


9.4 


16.0 


10 


-20.2 ± 


0.9 


2MASSW 


J2224438-015852 


L4.5 


18.02 


25.4 


22.2 


27.9 


3 


24.7 


20.9 


29.2 


6 


-37.4 ± 


3.4 


2MASSW 


J1507476-162738 


L5.0 


16.65 


27.1 


22.6 


31.6 


4 


27.2 


23.4 


32.7 


4 


-39.3 ± 


1.5 


2MASSI 


J0825196+211552 


L7.5 


19.22 










16.9 


11.3 


21.4 


6 


20.5 ± 


2.0 



The projected rotational velocities show no trend with spec- 
tral type (Fig. |3, as is also seen in the results from Mohanty & 
Basri ( 2003 ) over this spectral type range. There is similarly no 
correlation with I magnitude. 

6. Limits on rotation periods 

The projected rotation velocity, vsin/, is a lower limit to the 
equatorial rotation velocity, v. Assuming that a star rotates as 
a rigid sphere, v = InR/T, where R is the equatorial radius of 
the star and T is its rotation period. Given R, we can therefore 
derive the maximum rotation period, T max . More usefully, we 
can also derive the expected period and 'likely' range of pe- 
riods consistent with v sin i. This may then be compared with 
published periods obtained from photometric monitoring. 

Rigid body rotation is conventionally supported by the the- 
oretical argument that these cool stars are probably fully con- 
vective (e.g. Chabrier & Baraffe 1997 1 although counter argu- 
ments exist (Mullan & MacDonald 2001 1. Evolutionary mod- 
els show that a few hundred million years after formation, 
UCDs with masses between 0.04 M G and 0.09 M G all have radii 
around 0.1 R (Chabrier et al. Q997). Specifically, between 
0.5 Gyr and 1.0 Gyr the radii are in the range 0.10-0. 12 R Q , 
shrinking to 0.085-0.1 15 R at 5 Gyr. As the observed field 
UCDs reported here are probably at least 1 Gyr old (Gizis et al. 
2000 1, I adopt a radius of 0. 1 R for all targets when deriving 
period ranges. The effect of deviations from this are discussed 




2 4 

spectral type (0 = L0, 1 : 



= L1 
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etc.) 



Fig. 2. Projected rotational velocities for ultra cool dwarfs ob- 
served in this survey. The values are for the cross correla- 
tion against the L dwarf template, 2M1439+1929 (columns 
8, 9 and 10 of Table except for 2M1439+1929 itself and 
DEI 159+0057, for which results against the M dwarf tem- 
plates are shown. The cross symbol shows the expected value 
and the 'error bars' are the maximum and minimum values. 
Objects with equal spectral types are slightly offset on the spec- 
tral type axis for clarity. 

below, as are limits on R which can be derived from the data. 
In useful units, the relation between v sin /, T and R is 
R 1 



T/hrs =121.47 



0.1 R v/kms" 



(2) 
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Fig. 3. Derived period limits for the ultra cool dwarfs from the 
data given in Table [3]for the 2M1439+1929 template (except 
for 2M1439+1929 and DEI 159+0057, which are for the M 
dwarf templates). The star symbol is the expected period; the 
'error bars' show the maximum period and 90% confidence 
'minimum' period. 

Assuming that the spin axes of L dwarfs are randomly ori- 
ented in three-dimensional space, it can be seen that the prob- 
ability, P(/) di, of observing a spin axis between i and i + di is 
sin i di. From this we may simply derive that < sin i >- n/4 = 
0.785, and hence that < v > = 1.27 < vsin i > from which we 
calculate < T > from the above equation. The largest possible 
value of T, T max , occurs when sin / = 1.0, i.e. when the star is 
viewed edge on. The minimum value of T is arbitrarily small 
(as in the limit as sin i tends toward zero v tends toward infin- 
ity). A more useful statement concerns the probability, /, that 
v is not above some value, Vf. By fixing / at 0.9 we may then 
say that v is less than vt with 90% confidence. Generally, 

v/=(l-/ 2 r 1/2 vsini (3) 

so for / = 0.9, vj — 2.29 vsin i. The corresponding 90% confi- 
dence lower limit on the period, r m j n 9o, can then be determined 
from equation|2] 

These period limits are calculated from the vsin / values 
in Table [2] to give the period ranges reported in Table [3] and 
Fig -EJ the maximum period, r max , is calculated from the mini- 
mum value of v sin /, and the minimum period, TW^o, is calcu- 
lated using the maximum value of v sin i. The expected period, 
< T >, is calculated from the mean value of vsin/. Formally, 
the range 7 , m i n 9 U :7' max is the 90% confidence interval for the 
period, and is one-sided because T max is an absolute maxi- 
mum (for fixed radius). However, r m ; n 9o is a rather conservative 
lower limit to the period because the statistical projection argu- 
ment has been applied to the upper limit of v sin i, 

7. Assessment of published variability and period 
data: individual objects 

Of the 16 UCDs with vsin / measured in this work, nine have 
variability or period data published by one or more authors: 
Bailer- Jones & Mundt (2001 hereafter BJM); Bailer- Jones & 



Lamm (2003); Clarke et al. j2(J(J2al l2002bl hereafter COT); 
Gelino et al. J20021 hereafter G02); Gelino ( 120021 . A further 
five UCDs have variability/period data for which v sin / values 
have been published by others. The data for these 14 UCDs are 
now compared. All variability data is in an I band unless stated 
otherwise. 

2M0746+2000. This was identified by Reid et al. ( I200T1 
as a near equal mass binary with a separation of 0.22" (2.7 AU) 
through HST imaging, with a magnitude difference of AI=0.6. 
There is no evidence for binarity from my spectral data which is 
presumably dominated by the brighter component. COT iden- 
tify significant variability in their light curve data with an am- 
plitude of 0.007 mag and derive a period of about 3hrs, al- 
though they underline that this interpretation may be compli- 
cated by the binarity. It is also based on just a single dim- 
ming in the data (the source was only monitored for 6.5 hrs). 
Nonetheless, this period is consistent with the period limits in 
Table [3] G02 report a period of 3 1 hrs for this source with a 
confidence of almost 5cr. This is inconsistent with the much 
smaller maximum period of 5.7 hrs in Tabled This maximum 
is based on a cross correlation peak clearly widened above the 
the minimum v sin / threshold and there is good consistency be- 
tween the individual orders and the two templates. The 31 hour 
period of G02 can therefore be ruled out as the rotation pe- 
riod. The rotation period and v sin / could only be made con- 
sistent if the star had a radius of at least 0.6 R . This would 
require the star to be either much more massive than an L 
dwarf or very young, perhaps only a few million years from 
its birth line (Chabrier & Baraffe [19971 see also Joergens et 
al. [2003 1. However, in a detailed study based on kinematics, 
spectroscopy and photometry, Reid et al. ( 2000 1 conclude that 
2M0746+2000 has a mass between 0.07 M and 0.09 M and 
an age of at least 1 Gyr. A probable explanation for the 3 1 hour 
period from G02 is that it is not a real period (a 5cr detection 
is not large for a periodogram). Alternatively, it could be an 
alias of a shorter period. A third possibility is that this period 
is not related to the rotation at all, and may not even be a sta- 
ble period. It could instead be an artifact of unstable light curve 
modulation due to non-static surface features, as was suggested 
for some UCDs in BJM. Gelino ( 2002 1 also reports J,H,K mon- 
itoring observations of this star. Although there is formally no 
variability across the full light curves, there is strong evidence 
for a rise of 0.02 mag in J and K lasting 1.5 hrs (over a total 
time span of 5.5 hrs: see Fig. 4.5 in that work). 4 Interestingly, 
this rise is similar in duration to the dip seen in the same object 
by COT. 

2M1412+1633. G02 report a non-detection of variability 
on scales above 0.025 mag (1.4 times the reported RMS) on 
time scales of up to 85 days (although most data were obtained 
over 7 days according to Fig. 3.7 of Gelino 2002 1. Likely pe- 
riods from Table [5] could in principle have been detected, al- 
though the time sampling of G02 is very sparse so sensitivity 



4 Gelino states that this is not significant, but this conclusion ap- 
pears to be based on comparing the amplitude with the photometric 
error in a single measurement. Yet the 1.5 hrs rise is observed in about 
30 points systematically displaced from the mean, and is clearly a sig- 
nificant deviation. 
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Table 3. Expected, < T >, maximum, T max , and 'minimum', r m j n 9o, 
rotation periods for the ultra cool dwarfs taking into account both the 
statistical distribution of i and the range of v sin i given in Table. |2 
< T > is calculated using the mean v sin i and the expected value 
of sin i (= 7r/4). T max is derived from the minimum v sin i and max- 
imum sin i (= 1). r m i n 9o is calculated from the maximum v sin i and 
from that value of i (=26°) for which there is a 90% chance that 
the inclination lies above this: i.e., there is only a 10% chance that 
the period is shorter than r m j n 9o when adopting the maximum v sin i 
(and even less chance when adopting a more likely v sin z). 



Name 



rotation period / hrs 



M dwarf template 



< T > 



'- min90 



2M1439+1929 template 

< T > T mm on 7* max 



DE1431-1953 

DEI 159+0057 

2M0746+2000 

2M1412+1633 

2M1439+1929 

DEI 44 1-0945 

2M1145+2317 

2M1334+1940 

2M1029+1626 

DE1047-1815 

2M0913+1841 

SD1203+0015 

2M1615+3559 

2M2224-0158 

2M1507-1627 

2M0825+2115 



2.02 
2.72 
3.50 
5.53 
8.54 
5.50 
7.65 
3.80 
3.38 
5.87 
6.38 
3.47 
7.90 
3.77 
3.53 



1.01 
1.34 
1.73 
2.76 
4.14 
2.57 
3.68 
1.89 
1.63 
2.90 
3.01 
1.64 
3.84 
1.90 
1.68 



2.75 
3.81 
4.74 
8.26 

12.65 
8.21 

11.57 
5.50 
5.15 
8.10 

10.21 
5.04 

11.35 
5.47 
5.37 



2.58 
1.28 
3.71 
5.83 

6.02 
7.53 
3.77 
3.42 
6.38 
4.71 
3.02 
7.47 
3.87 
3.52 
5.66 



1.28 
0.64 
1.84 
2.76 

2.75 
3.58 
1.77 
1.60 
3.01 
2.35 
1.44 
3.32 
1.82 
1.62 
2.48 



3.60 
1.76 
5.28 
8.74 

9.13 
11.57 
5.65 
5.33 
9.96 
6.75 
4.74 
12.92 
5.81 
5.19 
10.75 



to the likely periods is low. If there were sufficient temporal 
sensitivity, then any net rotational modulation must have been 
smaller than 0.025 mag. 

2M1439+1929. This was not detected as variable by any 
of three independent monitoring programs: BJM derive an up- 
per limit of 0.01 mag over timescales of 1 hrs to 100 hrs; a 
slightly higher upper limit over longer timescales comes from 
G02. Bailer-Jones & Lamm (2003) restrict the J and K band 
variability to be less than 0.04 mags on time scales between 20 
min and 13 days. All of these programs had sampling dense 
enough to detect likely rotation periods from Table [3] I con- 
clude that the amplitude of rotational modulation is less than 
0.01 mag in I and 0.04 mag in J and K. 

2M1145+2317. While BJM reported a period of 
1 1 .2 + 0.8 hrs, a second monitoring epoch reported in the same 
paper confirmed the variability but not the period. This lack 
of stability led BJM to suggest that the modulating pattern on 
the star was not stable: if surface features were evolving dur- 
ing the second epoch on a timescale less than the rotation pe- 
riod, this could have 'masked' the rotation period. The present 
v sin i data suggest a period between 3.6 hrs and 1 1 .6 hrs, con- 
sistent with the 1 1.2 hrs period and this masking hypothesis. (It 
is not possible that the modulating features were simply weaker 
in the second epoch as then no variability would have been 
found.) A period of 11.2 + 0.8hrs corresponds to an equato- 
rial rotation speed of 10.8 ± 0.8 km/s (eqn.|2j. Combining this 
with the vsinz range of 10.5-14.8 km/s implies that the in- 



clination angle, i, lies in the range 62°-90°, with this lower 
limit being rather conservative. The expected value is z'=76°. 
Thus if the rotation period of 11.2 ± 0.8 hrs from BJM is cor- 
rect, 2M 1145+2317 would appear to be a near edge-on ro- 
tator (which makes detection of rotation by monitoring more 
likely). If the stellar radius were larger, this would decrease the 
expected and minimum inclination angles (and eventually the 
maximum too), and vice versa. Alternatively, we can use this 
period plus the minimum value of vsinz (10.5 km/s) to place 
a lower limit on the radius (i.e. corresponding to z'=90°) of this 
L1.5 dwarf. This gives R min = 0.097+0.007 R Q , the uncertainty 
arising from the uncertainty in the period. This is consistent 
with structure models for a large range of ages (see section|6j. 

2M1334+1940. The variability detection of this object 
was one of the most significant in BJM; the detected period 
of 2.68 ± 0.13 hrs was relatively significant (12<x). This is con- 
sistent with the present period deductions and corresponds to 
an equatorial rotation speed of 45.3 ± 2.2 km/s. Combining this 
with the vsinz range from Table|2]of 2 1.5-29. 9 km/s implies a 
range of z'=27°— 44°, with z'=34° the expected value. Although 
this star is quite faint, we can have some confidence in this 
given the agreement in v sin i values for the two templates and 
the reasonable significance of the period detection. As in the 
previous case, we can instead use the data to derive the min- 
imum radius, which is R m i n = 0.047 ± 0.002 R . This is a 
smaller limit than in the previous case, because these data per- 
mit smaller inclination angles for 2M1334+1940. 



C. A.L. Bailer- Jones: Rotation velocities of L dwarfs 



9 



2M1029+1626. G02 report no variability in this L2.5 
dwarf, although as this is based on only a few observations per 
night (Fig. 3.4 of Gelino 2002 1, this program would probably 
not have been sensitive to the short periods deduced in Tabled 

2M0913+1841. This was found to be variable by BJM 
but with no significant period. These observations should have 
been sensitive to the likely periods reported here (although 
plausibly the sampling in BJM could have thwarted its detec- 
tion: see Fig. 1 in that article). Assuming that this variability 
is due to changes in photospheric features, then these features 
must be changing in pattern and/or brightness at least as fast 
as the likely rotation periods (3-10 hrs) in order to mask a ro- 
tationally modulated signature. This is the conclusion of BJM 
based on the assumption that this L dwarf had a rotation period 
in the range of detectability. That this is now demonstrated to 
be likely adds some strength to that conclusion. 

SD1203+0015. As with 2M0913+1841, this was found to 
be a non-periodic variable by BJM. On the other hand, it was 
intensively monitored by Bailer- Jones & Lamm (2003 1 in the 
J and K bands, who found no variability above the larger limit 
of 0.04 mags over timescales of 20min to 7 days. The same 
comments as for 2M0913+1841 therefore apply. 

2M2224-0158. G02 report no variability, but as with 
2M1029+1626, there were only a few observations over each 
of a few nights, so periods of a few hours could have remained 
undetected. 

There are five other UCDs with published vsin/ values 
which also have variability data in the literature. These are now 
briefly discussed to complete this survey of UCD rotation and 
variability data. 

2M1146+2230 (L3V). This is an optical binary with a sep- 
aration of 0.29" (7. 6 AU) and AI=0.3mag (Reid et al. l200ll . 
Mohanty & Basri (2003) determine v sin i — 32.5 + 0.2 km/s, 
presumably for the brighter component. BJM detected 
marginally significant variability, with possible rotation peri- 
ods at 5.1±0.1hrs (15(f) and 3.0 + 0.1hrs (6(f). The former 
is slightly too long for the measured v sin i, although it would 
be consistent if the radius were 0.1 3R Q (more if the star is not 
observed edge on). This radius is consonant with the evolu- 
tionary models, especially if this L dwarf is less than 1 Gyr in 
age (Chabrier & Baraffe fT997> . The latter period would also be 
consistent, but is of too marginal significance to warrant further 
discussion. Neither G02 nor COT detected variability in this 
source, although the upper limits of both are consistent with 
the BJM detection. Given the sampling, either BJM or COT 
could have detected the likely periods, but the weight of evi- 
dence is that the modulation contrast was not large enough (or 
not stable enough) to permit detection. 

LP944-20 (M9V). vsin/ was determined to be 
28.3 + 2 km/s by Tinney & Reid (19983 and 39 + 2 km/s 
by Mohanty & Basri d^DOll . Tinney & Tolley 0999) detected 
variability with an amplitude of 0.04 mag and timescale of 
2 hrs. (With a p value of only 0.02 this would not have been 
counted as variable on the criteria of BJM or G02, although 
the fact that the detection is made in two narrow bands 
compensates for this.) Although these authors did not associate 
this with a rotation period, it is consistent with the maximum 



rotation periods of 4.3 hrs and 3.1 hrs implied by the above 
vsin/ values. 

DE1228-1547 (L4.5V). Tinney & Tolley ( fT999b did not 
detect variability in this, although their observations were in 
principle sensitive to the periods commensurable with the v sin i 
of 22 km/s reported by Basri et al. ( 2000 ) and Mohanty & Basri 
(2003). This implies low contrast surface features. 

BRI0021-0214 (M9.5V). Tinney & Reid ( fT998l estab- 
lished vsin/ = 42± 8 km/s for this compared to 34 ±2 km/s 
from Mohanty & Basri (2003). Martin et al. (1200 lb detect sev- 
eral periods in this source from I band monitoring, in particular 
at 4.8 hrs and 20 hrs. Both are too long for the measured v sin i 
range unless the radius is at least 0.13 R and 0.5 R Q respec- 
tively - more if the star is not observed edge on - something 
which Martin et al. rule out. This is taken as evidence for the 
masking hypothesis. 

Kelu-1 (L2V). This is a very rapid rotator, as reported 
both by Basri et al. ( 120001 v sin / = 60 ± 5 km/s) and Mohanty 
& Basri (2003J vsin/ = 60 + 2 km/s). Clarke et al. (2002a) de- 
tected Kelu-1 to be periodically variable with a period of 
1 .80 ± 0.05 hrs and peak-to-peak magnitude of 0.01 1 mag. This 
is consistent with the v sin / values, and is discussed at length 
by Clarke et al. 

8. Conclusions 

I have presented vsin / measurements for 16 ultra cool field 
dwarfs with spectral types between M9 and L7.5. 14 of these 
are new, doubling the number of L dwarfs with measured v sin i. 
14 UCDs have vsin / between 10 and 30km/s, the other two 
between 30 and 40 km/s. This confirms a previously published 
result that, compared to the M dwarfs, L dwarfs are fast rota- 
tors, with implications for their angular momentum evolution. 
In addition to establishing mean v sin / values I determined a 
conservative v sin i range for each UCD which encompasses the 
asymmetric uncertainties in fitting the cross correlation peak: 
this full range is typically 7 km/s. 

With theoretically predicted radii of 0. 1 R , the expected 
rotation periods of these UCDs (i.e. using < sin / >= 7r/4) is 
3-10 hrs. The measured vsin / determines the maximum rota- 
tion period. Using the statistical distribution of /, a one-sided 
90% confidence range on the period of each UCD can be estab- 
lished: this is 5.3-12.1hrs and 1.8^1.0hrs for vsin/=10km/s 
and 30 km/s respectively. A larger radius for a given vsin/ 
would imply longer periods. 

From this analysis, suggested periods for UCDs from 
published monitoring programs were assessed. The follow- 
ing conclusions were drawn. The period of 11.2 + 0.8 hrs for 
2M 1145+2317 detected by Bailer-Jones & Mundt d2"0"0"fl is 
likely to be the rotation period for this LI. 5 dwarf, and the 
inclination angle is constrained to be /=62°-90°, i.e. a near 
edge-on rotator. Likewise, the period of 2.7 ±0.1 hrs detected 
for the LI. 5 dwarf 2M1334+1940 by the same authors is 
confirmed as the likely rotation period, and the inclination 
is /=27°-44°. The period of 31 hrs detected by Gelino et al. 
( 2002 ) for the brighter component of the binary 2M0746+2000 
can be strongly ruled out as a rotation period. In three cases 
(2M0913+1841, 2M1 145+2317 and SD1203+0015), variabil- 
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ity but no period was detected by monitoring programs, even 
though the present work derives likely periods which should 
have been detectable by the programs. Although all of these ob- 
jects are faint, this seems to support the 'masking hypothesis', 
the idea that surface features evolving in distribution and/or 
brightness faster than the rotation period could mask such a 
period from being detected by monitoring. The challenge to 
ongoing work remains to characterise this evolution. 
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